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Abstract 

We study a neutrino mass texture which can explain the neutrino oscillation data and 
also saturate the upper bound of the CP asymmetry Si in the leptogenesis. We consider 
the thermal and non-thermal leptogenesis based on the right-handed neutrino decay in 
this model. A lower bound of the reheating temperature required for the explanation of 
the baryon number asymmetry is estimated as O(10®)GeV for the thermal leptogenesis 
and O(10®)GeV for the non-thermal one. It can be lower than the upper bound of the 
reheating temperature imposed by the cosmological gravitino problem. An example of 
the construction of the discussed texture is also presented. 
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1 Introduction 


The discovery of the neutrino masses [1] gives large impact to the study of particle physics 
and astroparticle physics. In particular, it presents an interesting approach to the study 
of the origin of the baryon number (B) asymmetry in the universe, which is one of the 
important questions in these helds. The leptogenesis [2] based on the B — L violation due 
to the neutrino masses is considered to be the most promising scenario for the generation 
of the B asymmetry. During the recent few years, the leptogenesis based on the CP 
asymmetric decay of the heavy right-handed neutrinos [3] whose existence is required by 
the seesaw mechanism [4] has been extensively studied [5, 6]. 

Since the intermediate scale is generally necessary for the seesaw mechanism, it seems 
to be natural to consider the leptogenesis in the supersymmetric framework to guarantee 
the stability of that scale against the radiative correction. However, if we consider it in 
such a framework, a crucial problem called the cosmological gravitino problem is caused 
in relation to the generation of the right-handed neutrinos. If the reheating temperature 
Tr required to produce a sufficient amount of the right-handed neutrinos is a high value 
such as Tr > 10®GeV, the gravitino can be produced too much and its late time decay 
may disturb the nucleosynthesis [7].^ Thus, the production mechanism of the heavy right- 
handed neutrinos is the important ingredient for this problem. Several solutions for this 
difficulty have been proposed by now [9, 10]. 

On the other hand, the CP asymmetry [3] in the decay of the right-handed neutrinos 
is another crucial factor which plays the essential role to determine the generated lepton 
number (L) asymmetry. Its magnitude depends on the structure of the neutrino mass 
matrix, which is severely constrained to explain the neutrino oscillation data [1]. From 
a viewpoint of the leptogenesis, it is favorable that the neutrino mass matrix can realize 
the maximum value of the CP asymmetry [3, 11] or enhance its value [12, 13]. Thus it 
is important for the quantitative study of the leptogenesis to construct such a concrete 
model for the neutrino mass matrix as done in [5, 14, 15] and to proceed the investigation 
based on it. In this paper we present an example of the neutrino mass matrix and estimate 
the reheating temperature required to produce the sufficient B asymmetry. 

The paper is organized as follows. In section 2 we present a neutrino mass texture 
^If the gravitino is the lightest superparticle as in the gauge mediated supersymmetry breaking, there 
is no gravitino problem even in the case of Tr = O(10^°)GeV [8]. 
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and discuss its phenomenological features. An example for its construction is discussed 
in appendix A. In section 3 we apply this model to both the thermal and non-thermal 
leptogenesis. We discuss the lower bound of the reheating temperature required for the 
production of the sufficient B asymmetry. Section 4 is devoted to the summary. 


2 Neutrino mass texture 


We consider the minimal supersymmetric standard model (MSSM) extended with gauge 
singlet chiral superhelds W which correspond to three generation right-handed neutrinos. 
An effective superpotential for the neutrino sector is assumed as follows: 

W=j2 , ( 1 ) 

i,j=l ^ Z J 


where Li and H -2 are the lepton doublet and Higgs doublet chiral superhelds, respectively. 
In this paper we use the same notation for both a superheld and its component helds. 
The right-handed neutrino mass matrix M. and the Dirac mass matrix tud induced from 
the hrst term in this superpotential are assumed to take the form as^ 


M = 


Ml 

m 

0 


m 0 ^ 
M 2 0 
0 M 3 ! 


mo = h''{H2) 


^ 0 0 0 ^ 

a a' 0 

^0 b b' j 


( 2 ) 


where the charged lepton mass matrix is considered to be diagonal. Although each element 
of M. is supposed to be real, is assumed to be a complex matrix. 

If the hierarchical structure 


m. Ml M2 M3 ( 3 ) 

is assumed in the right-handed neutrino sector, the eigenvalues Mj of the mass matrix M. 
can be approximated to 

(a) Mi(~Mi), M 2 , M 3 , (form 2 <MiM 2 ), 

(b) Mi(~ M 2 sin 2 e), M 2 , M 3 , (for m2 > Ml M 2 ), 

^This model can be considered as a simple extension of the model in [16]. It is realized by adding a 
right-handed neutrino Ni to the original one in such a way that Ni weakly couples to N 2 alone. 
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where we use Mi ~ Mi {i = 2,3) and sin,^ ~ m/M2. These two cases are studied in 
the following part. The structure of Ai and mn can be effectively realized by imposing 
a snitable symmetry on the snperpotential at the high energy scales. We give such an 
example for the construction of Ai and mr, in appendix A. 

If we change W into the Ai diagonal basis W, the Dirac neutrino mass matrix is 
transformed into 


rfiD = h''{H2) 


( 


—asin,^ 
a cos^ 


—a' sin ^ 
a' cos ^ 


0 ^ 
0 


(5) 


V 


J 


0 b b' 

Applying the seesaw mechanism to these matrices, we can obtain the light nentrino mass 
eigenvalues and the MNS matrix. Here, for the simplicity, we put a' = y/2a and U = b,^ 
and then the light neutrino mass eigenvalnes are fonnd to be 


mi = 0 , m2 — 


2\a\' 

Wo 


mo, ~ 


m 

Mo 


( 6 ) 


The MNS matrix has the bi-large mixing form snch as 


( T2 


[/ ~ 


V 


—i icos6* ^cos6* 
I —^cos0 ^cos6* 


(7) 


where we neglect the CP phase in this expression. 

Now we can compare these results with the present experimental data. Since the 
neutrino mass eigenvalnes are assumed to be hierarchical, the analysis for the neutrino 
oscillation experiments reqnires [1] 


2|a| 


- ^ (7 X 10 ^ ~ V Ami^ ~ (2 x 10 ^ (8) 


M2 

and sin 9 should satisfy 


M- 


sin 9 ~ 




|a|VM 2 1 _ 

v^|6|VM3 “72^ Amh^ 


0 . 1 . 


(9) 


This is consistent with the constraint sin0 < 0.16 which is imposed by the CHOOZ 
experiment [17]. The effective mass for the nentrinoless donble /3-decay is estimated in 

^We adopt these relations in the following study. Under this assumption the model contains seven real 
parameters. We can loose these strict equalities without changing largely the qualitative results given in 
this paper. 
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this model as 


mee 


< 


1 

2 



~ 2 X 10“®. 


( 10 ) 


It seems to be difficult to reach such a value in the next generation experiment. 

The lepton flavor violating processes such as p —67 can constrain the model. It 
has been suggested that these processes could impose the strong constraint because of 
the renormalization effect on the soft supersymmetry (SUSY) breaking parameters due 
to the off-diagonal Yukawa couplings. It can be very severe even in the case of the 
universal SUSY breaking in the gravity mediation scenario [18]. Here in order to And 
the conservative condition, we consider the universal soft SUSY breaking in the gravity 
mediation. The branching ratio of the flavor changing process ii —>■ ij'y is estimated by 
taking account of the one-loop contribution as [18] 


Br{i, ^ £, 7 ) = 


taffi /5 


-1 


87r2 


(StUq -|- Al)(hGh'')ij In 


M 


X 


M 


( 11 ) 


where mo, Aq and m^ represent the soft scalar mass, the SUSY breaking A parameter and 
the relevant slepton mass, respectively. Mx stands for the unification scale and M is the 
right-handed neutrino mass scale. Since Mi is irrelevant to the light neutrino masses as 
shown in eq. ( 6 ), M is appropriate to be taken as M 2 in the present model. 

If we assume mo — m^ ~ Ho and use eq. (5) for the Yukawa couplings K', each 
branching ratio for p —> eq and r —pq is estimated as^ 

_oi Afn / Mx\‘^ 0 ^ _11 

Br(u —> eq) ~ 3 x 10 ^ j (In —— ) tan^ /? < 1.2 x 10 , 

mo V M 2 J 

Hr(r ^pq) ~4 X 10“=^°^ fln^) tan^/? < 1.1 x 10“®. (12) 

mo V M 2 J 

If we take Mx = 10^® GeV and mo = 100 GeV, for example, the present experimental 
bounds can be satisfied for M 2 ^ lO^^GeV and M 3 < 10^^ GeV even in the case of large 
tan/? such as tan/? ~ 50. This means that no lepton flavor violating decays ^i —> £^7 
contradict with the present model as far as the universal SUSY breaking is assumed even 
in the gravity mediation scenario. 

A remarkable feature of the model is that there are no constraints on Mi and sin,^ 
from the neutrino oscillation data and other present available experiments. If we apply 
■^The decay r ^ eq is automatically forbidden as a result of the present texture of the neutrino mass 
matrix. 
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this neutrino mass texture to the leptogenesis, Mi and sin ^ may be constrained to explain 
the B asymmetry. In the next section we focus our study on this point. 


3 Application to leptogenesis 


The decay of the heavy right-handed neutrinos can produce the B — L asymmetry. Then it 
may explain the B asymmetry in the universe since the sphaleron interaction can convert 
a part of the B — L asymmetry into the B asymmetry [19]. The L asymmetry or the 
B — L asymmetry induced through this decay is produced as a result of the CP asymmetry 
caused by the interference between the tree and one-loop diagrams. 

The CP asymmetry appeared in the iVj decay can be generally expressed as [3] 


£i 


El T{N, ^ l,H2) - El r(iV, ^ L,H2) 

El t { n , ^ L1M2) + El r(iv, ^ I1M2) 


1 1 

Stt 


a kj^i 



(13) 


where f{x) contains the contributions from both the vertex correction and the self-energy 
correction, and it has an expression 



+ 


2 ■ 

X — 1. ' 


(14) 


Applying this formula to our model in which the hierarchical structure for the right-handed 
neutrino masses is assumed, we obtain 


£1 ~ 


^3 ^ 




Stt sin f3 


sin 2 x. 


£2 - 


Am^trvMo 


Stt ^ 


■sin 2 x, 


1 /Ml M3 2 . , M2 , Ms 


IGtt sin^ (3 


ln^sin^/-h ^In^cos^/1 sin2x, (15) 

Ms Ml Ms M2 7 ^ ^ 


where (M 2 ) = nsin/? and y = arg(a* 6 ). The formulas in eq. (15) show that £3 can be 
much smaller than £ 12 . The interesting point of this result is that ei is almost equal to 
the expression which saturates its upper bound given in [11].® Thus the present texture 
for the neutrino masses seems to be favorable to induce the B asymmetry. 

model with this feature has been discussed in [14] already. However, it seems not to have a 
satisfactory structure for the explanation of the neutrino oscillation data. 
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By using these expression for Sj, the L asymmetry resulting from the decay of the 
right-handed neutrinos can be estimated. If we put the excess of the number density of 
the right-handed neutrino iVj from the equilibrium one as n* and the entropy density in 
the comoving volume at the latest Ni decay as s, the produced L asymmetry through this 
decay can be expressed as 


Yl = - ^ —SiKi, 




(16) 


where Ki represents the washout effect which depends on the strength of the Yukawa 
couplings in eq. (5). The sphaleron interaction which is in the thermal equilibrium at 
the temperature 10^ GeV ^ Tgph ^ 10^^ GeV converts a part of the B — L asymmetry into 
the B asymmetry in such a way as ub/s = —(8/15)(nL/s) in the MSSM case [19, 20]. Thus 
ul/s should satisfy |ni/s| ^ 10“^° to realize the observed value ub/s ~ (0.6 — 1) x 10“^°. 

Since the number density n* of the right-handed neutrino Ni depends on its generation 
mechanism, we need to £x it for the quantitative estimation of ul/s. In the following 
study we consider both the thermal and non-thermal scenarios for their generation. In 
the non-thermal leptogenesis we mainly consider that the right-handed neutrinos couple 
to the inflaton directly and then the inflaton decay produces the right-handed neutrinos 
non-thermally. 


3.1 Thermal leptogenesis 


Since the right-handed neutrino masses are supposed to be hierarchical in the present 
model, the L asymmetry is expected to be produced as a result of the out of equilibrium 
decay of the lightest right-handed neutrino Ni as studied in a lot of works [5, 6]. In fact, 
it can be easily checked in the present model. If we use eqs. (5) and (8), the decay width 
of Ni can be estimated as 


F - ~ _ 

167r 


AmLMiM2 


sin^ f3 

^ _3 

IbTT sin^ (3 


sin^^ 


Ml 


cos^ ^ ~ lO'" 


^107 GeVy 

/ M 2 
UQio GeV 


Mo 


IQio GeV 


x2 

J cos^^. 




F - ~ 

87r 


1 ^ / M3 

Itt v'^sin^P V10^^ GeV/ ’ 


(17) 


where these decay widths are given in the GeV unit. This shows that T ^N 2 ^ ^N 3 
is satished. The asymmetry generated by the decay of N 2,3 is washed out through the L 
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violating scattering mediated by the thermal Ni etc. and then ^ 2,3 ^ 1- 

If we use the thermal number density of the relativistic particle for rii and s = 
in eq.(16), the L asymmetry generated through the Ni decay is expressed as 


riL 1 

- ~ —EiKi, 

s g^f 


( 18 ) 


where g^ is a degree of freedom for the relativistic particles at this period and g^ ~ 200 
in the MSSM. In this expression we should note that Ki includes also the efficiency factor 
to generate Ni in the thermal bath other than the washout effect since we suppose that 
there are no thermal right-handed neutrinos initially. 

As is well known in the thermal leptogenesis [6], there is an important quantity mi 
which is related to ki and then the strength of the relevant Yukawa couplings of Ni. 
It controls how many Ni is produced in the thermal equilibrium and also how much L 
asymmetry is washed out. In the present model mi is expressed as 


fhi = 


11 


3 /T ~ ^2 . 2 f 


Ml 




for (a), 
for (b). 


(19) 


As is found from eqs. (15) and (19), ei and rhi can be independent from each other because 
of the freedom of sin,^. We may expect that there is generally some correlation between 
these parameters from their dehnitions (13) and (19). However, the special texture can 
make them independent in the present model. This feature may cause a substantial 
influence on the reheating temperature required for the leptogenesis. 

If we use the formula in eq. (15), the CP asymmetry ei required from the B asymmetry 
in the universe is estimated as 


kil ^ 10 


-8 


X 


Ml 


108 GeV 


10 " 


( 20 ) 


where g^ ~ 200 is used. In this estimation the maximum CP phase |sin2x| ~ 1 and 
sin j3 = 1 are also assumed.® From this condition we obtain Mi ^ 10® GeV for the lower 
bound of the Ni mass, which is the ordinary result in the out of equilibrium decay of the 
thermally produced Ni. On the other hand, the effective mass fhi is estimated as 


mi 

10-2 eV “ 


^2 . 2 f r ( \ 

— sm ^ for (a), 

1 for (b). 


( 21 ) 


®This assumption for sin j3 brings no crucial difference as far as tan (3 is in the interesting region such 


as 1 < tan/3 < 50. 



While fhi takes a fixed value in the case (b), there is the freedom sin,^ to tune rfii into 
the desirable value in the case (a). Eq. (21) shows that the efficiency factor can be in 
the favorable region for the leptogenesis in the case (a) but it seems to be larger than the 
favorable value in the case (b) [5, 6 , 15]. 

The necessary condition for the out of equilibrium decay of Ni is given by > Tj^^. If 
we use eq.(17) for this condition, the successful leptogenesis requires that the temperature 
T at the period of the Ni decay should satisfy 

T > Tmin — \lM 1 M 2 Sin ( 22 ) 

In the case (b) we find T^ain — Mi and then Tr > Mi should be satisfied. Thus we expect 
the similar result for the efficiency factor to the previous works [5, 6 ]. On the other hand, 
since the case (a) is realized for Mi > M 2 sin^ we find that T < Mi could be consistent 
with the condition for the out of equilibrium decay. Such a situation seems to be realized 
for a sufficiently small sin,^ without conflicting with the neutrino oscillation data. The 
small sin,^ can also make the washout effect negligible. Although this seems to suggest 
the possibility that the reheating temperature Tr may not be necessary to be high enough 
compared with Mi, however, the sufficient number of Ni may not be produced due to the 
Boltzmann suppression in the thermal equilibrium. We need to solve a set of Boltzmann 
equations numerically for the quantitative study of the relation between Tr and Mi. 

We study these points by solving numerically a set of Boltzmann equations for the 
MSSM presented in [5]. If we use eq. (8) and assume | sin 2x| = 1, the model parameters in 
this calculation are Mi 2,3 and sin,^. As an initial condition for the Boltzmann equations, 
we assume that both the number density of W and the L asymmetry are zero at zq = 0. 
A dimensionless parameter is defined as z = Mi/T. The temperature corresponding to 
Zq may be considered to correspond to the reheating temperature Tr. 

In Fig. 1 we give a solution of the Boltzmann equations with zo = 0.01 and the input 
parameters such as Mi = lO^GeV, M 2 = 10^°GeV, M 3 = lO^^GeV and sin,^ = 0.02. 
This corresponds to the case (a) with Tr = lO^^GeV. The case (b) cannot yield the 
sufficient L asymmetry. This can be understood as follows. In this case sin,^ is required 
to take a rather large value to realize Mi < M 2 sin^ ^ satisfying both constraints such 
as Ml > 10®GeV and M 3 < lO^^GeV, which are imposed by the previously discussed 
phenomenological constraints. Such a sin,^ makes the Ni Yukawa couplings larger and 
then the washout effect becomes effective. This is also suggested by eq. (21). 
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Fig. 1 A typical solution of the Boltzmann equations in the case of the thermal generation of Ni. We 
define Yi as Yi = rii/s (i = Ni,N2,L) and Y±i = {nsN^ ihsArJ/s where SNi stands for the sneutrinos. 
Y^^ is the value in the equilibrium. 

In Fig. 2 we show the L asymmetry \Yl\ as a function of sin,^. In the left panel we plot 
\Yl\ for the various values of Mi. In the right panel \Yi\ is plotted for the various values of 
Zq. If we use eqs. (21) and (22), we hud that the input parameters adopted to draw Fig. 2 
give ?fii ~ 10“^~“^eV and Tr > lO^GeV. By using this kind of hgures, we can search the 
lower bounds of Mi and Tr required to explain the observed B asymmetry for the hxed 
^ 2 , 3 - We practice this analysis changing the values of M 2^3 within the allowed region 
discussed in the previous part. As the result, for the explanation of the B asymmetry 
based on the present model, we hud that the lower bounds of Mi and Tr can be estimated 
as 

Ml > 5 X 10® GeV, Tr > 6 x 10® GeV, (23) 

and also sin,^ should be 0(10“^). 

Although the obtained lower bound of the reheating temperature is comparable with 
the lowest value discussed in other neutrino mass matrix models where Tr ~ 10®“^°GeV 
is usually suggested, we cannot make it much lower. Since the lower bound of the mass 
eigenvalue Mi is determined by eq. (20), this result seems not to be avoided in the 
thermal leptogenesis as far as we do not assume the degeneracy among the masses of the 
right-handed neutrinos. Recently, in [21] the upper bound for the reheating temperature 
required from the gravitino problem is estimated as 10^“^GeV if the gravitino has the 
mass in the range 10^“®GeV. If we do not suppose the light gravitino scenario and we 
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0 0.05 0.1 0 0.05 0.1 

sin ^ sin ^ 


Fig. 2 The L asymmetry |Yl| as a function of sin^. Horizontal thin lines represent the desirable region 
to explain the observed B asymmetry. In the left panel Mi is varied keeping other parameters fixed in 
such a way as M 2 = lO^^GeV, M 3 = lO^^GeV and zq = 0.01. In the right panel we vary the zq value 
keeping others fixed as Mi = 10®GeV, M 2 = lO^^GeV and M 3 = lO^^GeV. 

follow this bound, the present model is unable to be reconciled with the gravitino problem. 
We need to consider the initial Ni to be yielded in other way. As such a possibility, we 
study the non-thermal leptogenesis in the next part. 


3.2 Nonthermal leptogenesis 


In this subsection we consider that the right-handed neutrinos W are produced through 
the decay of the inflaton. This kind of model has been discussed in [9]. The interaction 
between the inflaton superheld and W is assumed to be given by the snperpotential 

3 

W = ^ (24) 

i=l 

After the inflation ends, the inflaton 0 starts to oscillate and decays to reheat the nniverse 
into the temperatnre Tr. A part of its oscillation energy p of the inflaton is converted into 
Ni throngh its decay a.t H czl V The decay width T,^ of the inflaton 0 can be expressed 
as 

n = + (25) 

where the ellipses stand for the contribution from other decay modes of the inflaton and 
we assnme them to be negligible. The coupling constants A* are constrained in snch a 
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which is derived from the condition ~ F^. From this we hnd that the couplings between 
the right-handed neutrinos and the inflaton can be small enough not to affect the inflaton 
potential. The inflaton mass can depend on the assumed inflation model. However, it 
should satisfy > Mj to guarantee the inflaton decay into the right-handed neutrinos 

N^. 

If we use Bi to denote the branching ratio for the decay 0 — N^, we have the energy 
relation pBi = MiUi where p = Thus the non-thermally generated number 

density n* of iVj can be written as 

^ _ STuBj 

s AMi ' 

As we mentioned below eq. (15), £3 is much smaller than £ 1 ^ 2 - In the case dominated by B^ 
the produced L asymmetry is expected to be very small unless Tr is rather large. Then we 
concentrate our study into the case dominated by Bi^ 2 - Since we hnd ei/Mi = € 2 /M 2 from 
eq. (15), the L asymmetry generated through the immediate decay of Ni ^2 is estimated 
by using eq. (16) as 


nL 

s 


^ , .T/ (/s:ii?i + /.2i?2)sin^2x 

IbTT sm p 

~ ( lO«^Gev ) 


(28) 


Eq. (28) suggests that Tr should be larger than 10® GeV to explain the B asymmetry 
in any case. Moreover, we can hnd that the Mi dependence of rii/s is conhned into the 
washout factor Ki. If we consider that Mj is larger than Tr, the washout ehect due to iVj 
is expected to be suppressed by the Boltzmann factor. Thus for the larger Mj in such a 
region, ul/s becomes insensitive for the change of the values of Mi. 

The estimation of the L asymmetry in eq. (28) is justihed only if Ni decays into the 
light helds immediately after its production [9]. This requires that FT ~ T,^ < T^y, should 
be satished for all iVj which have the substantial branching ratio Since the inflaton 
decay width can be estimated by using iF ~ T^ as 

r* - (29) 

tWpl 

^In the construction of the mass matrices presented in appendix A, we find that ^ i?i is satisfied 
if the inflaton has no global charges. However, we consider the general case here. 
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Fig. 3 A typical solution of the Boltzmann equations in the case of the non-thermal generation of Ni. 
The definitions of Yi are the same as the ones in Fig. 1. 


we can write the condition for the justihcation of eq. (28) by applying eqs. (26) and (29) 
to ^ in the form as 


Tr y 
106 Gev; 


106 


Ml \ 

107 GeV j 


Ms 

1010 GeV 



(30) 


Since this condition can be easily satished for the desirable values of Tr, Mi and Ms, we 
hnd that eq. (28) can be validated in our interested case. However, it is also possible that 
the immediate decay condition is not satished for Ni. This occurs for Tj^^ < T^ < T^^ in 
the case of Bi ~ i?s- In that case we should take account that the L asymmetry produced 
through the iVs decay may be washed out by the late entropy release due to the Ni decay 
other than by the usual thermal washout. This effect is discussed in appendix B. We 
should also check that the condition (30) can be consistent with the above mentioned 
condition > M^,. This consistency can be easily checked by applying eq. (26) to T^. 

In order to study the relation between Tr and M*, it is useful to classify the situation 
into three cases: (i) Tr < Mi, (ii) Mi ^ Tr ^ M 2 , (iii) M 2 ^ Tr. Among these three 
cases, both the cases (i) and (ii) can easily satisfy the condition (30). On the other hand, 
the case (iii) satishes it only for M 2 ^ lOMi sin^.^, which requires sin,^ 3> 0.1. Thus only 
the two cases (i) and (ii) seem to be promising for the leptogenesis consistent with the 
low reheating temperature. In fact, Tr ~ IO 6 GeV may be allowed in these cases with 
-Bi ,2 = 0 ( 1 ) and | sin 2 x| = 0 ( 1 ) if Ri or k ,2 can be 0 ( 1 ). 

The washout effect is expected to be mainly caused by the L violating interactions 
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due to the thermal iVj. Since there is the Boltzmann suppression for these processes in 
the case (i), their substantial washout effect cannot be expected. On the other hand, 
in the case (ii) Ni can contribute to the washout of the L asymmetry since there is no 
large Boltzmann suppression. To escape this situation the smallness of the Ni Yukawa 
couplings is required. This may be realized for the small sin,^ case. 

To take account of the washout effect quantitatively, we need to solve the Boltzmann 
eqnations nnmerically by using eq. (27) as the initial valne for the n* at Zq = Mi/Tr. 
In Fig. 3 we show a typical solution for the Boltzmann eqnations in the case (a). In 
this figure we assume | sin2x| = 1, Bi = B 2 = 0.5 and Tr = 3 x lO^GeV. The input 
parameters are taken as Mi = lO'^GeV, M 2 = lO^GeV, M 3 = lO^^GeV and sin,^ = 0.01. 
This hgure shows that the nnmber density of N 2 rapidly decreases following the Boltzmann 
distribntion. The L asymmetry reaches the hnal valne faster compared with the thermal 
case. In the case (b) the snfficient amount of the L asymmetry cannot be prodnced. The 
reason is considered to be the same as the thermal case. 

In Fig. 4 we plot the L asymmetry |Yr| as the fnnction of sin,^ for various values of the 
input parameters. We calculate Yl for the typical three models with different branching 
ratios and plot them with the different symbols, that is, the squares for i?i = 0 , i ?2 = 1 , the 
circles for Bi = B 2 = 0.5 and the triangles for Bi = 1, i ?2 = 0. The reheating temperatnre 
is assumed to be Tr = 3 x 10®GeV. The left panel corresponds to the case (i). In this 
hgure, as the input parameters we use Mi = 10®GeV, M 2 = lO^GeV, M 3 = lO^^GeV for 
three types of the black symbols and Mi = 10®GeV, M 2 = 5 x 10®GeV, M 3 = lO^^GeV 
for the white symbols. The right panel corresponds to the case (ii). In this case, as the 
inpnt parameters we use Mi = lO^GeV, M 2 = 5 x 10®GeV, M 3 = lO^^GeV for the black 
symbols and Mi = lO^GeV, M 2 = 10®GeV, M 3 = lO^^GeV for the white symbols. The 
typical featnre in these cases is that the sin,^ value can be smaller compared with the 
thermal case since we need not to prodnce W thermally. 

In both panels of Fig. 4 the larger \Yl\ is realized for the larger B 2 since the washont 
effect is smaller compared with the smaller B 2 case. If we make M 2 larger keeping Mi 
hxed in these hgnres, \Yi\ becomes a little bit larger bnt it seems to reach almost the npper 
bonnd in this setting. In the left panel the condition (30) is satished only for sin,^ ^ 10“^. 
In the case of i ?2 = 1, however, this condition shonld be replaced by F,^ < F^^ and it is 
satisfied for all region of sin,^ in this hgure. In the case of Bi = B 2 = 0.5, for sin,^ < 10“^ 
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sin ^ sin ^ 

Fig. 4 The L asymmetry Yl as a function of sin^. Horizontal thin lines represent the desirable region 
to explain the observed B asymmetry. The explanations for the symbols are presented in the text. 


we should take account of the additional washout effect to the L asymmetry produced 
through the N 2 decay, which is discussed in appendix B. It introduces the suppression 
factor \l Ml M 2 sin,^/TR. It takes a value smaller than 0(10“^) in the present case. On 
the other hand, the L asymmetry produced by the late Ni decay cannot be a sufficient 
amount because of the low reheating temperature. Thus we find that the models with the 
substantial Bi cannot explain the B asymmetry for sin,^ < 10“^. The Yukawa couplings of 
Ni and N 2 are proportional to sin,^ and cos^, respectively. This fact affects the behavior 
of \Yi\ as the function sin,^. In fact, the figures show the slight increase in the case of 
B 2 = 1 and the decrease in the case of i?i = 1 when the sin,^ value becomes larger. 

In the right panel the condition (30) is satisfied only for sin,^ > 0.1. The situation 
is the same as the left panel in the case B 2 = 1. Thus in the case (ii) only the model 
with i ?2 ~ 1 can have the possibility to explain the B asymmetry. The magnitude of \Yi\ 
becomes smaller for the larger sin,^ even in the case of B 2 = 1. Since Ni can be produced 
thermally, the large sin,^ makes the washout effect more effective. This feature explains 
the \Yl\ behavior against sin,^ in the right panel. 

We practice this kind of analysis changing the input parameters within the allowed 
region. As a result of such study, we find that the lower bound of the required reheating 
temperature to explain the B asymmetry is Tr ~ 3 x 10® GeV for both cases (i) and (ii). 

Finally, we briefly comment on other possibility for the non-thermal leptogenesis [10]. 
In the early universe the scalar potential of the sneutrino Ni may be flat enough to deviate 
largely from its potential minimum.® If this happens during the inflation, the condensate 
®The sneutrino Ni can be an inflaton itself as discussed in [22]. However, we do not consider this 
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of Ni starts to oscillate at ~ Mi and decays at M ~ This oscillation may 

dominate the energy density of the universe at a certain time after the reheating due to 
the inflaton decay because of its behavior as a matterWe assume that it is the case 
here. 

Since its energy density is expressed by = M^|iVip, the Ni number density ni is 
estimated as Mi|Wp. Thus the ratio of the L asymmetry produced through its decay to 
the entropy density is estimated as 


s 


2Mi|iVi|2 _3Tr 

- 

s 2 Ml 


(31) 


where in the last equality we use the above mentioned assumption Pf;^ = for the 

energy density. If we use eq. (15) for the CP asymmetry ei, we obtain the similar result 
for ul/s to the previous example. However, in this case T^ > T^^ should be satished 
and this condition imposes Tr > \jlOM 1 M 2 sin,^. Thus the expected L asymmetry is 
estimated as _ 


^ ~ 10 - 1 ° 


T, 


R 


10 


\/MiM 2 sin ^ 


s 106 GeV ~ 105 GeV ’ 

where we assume |sin2x| = 1 and Ki = 1. This relation gives a constraint for the 
undetermined parameters in the present neutrino mass texture based on the B asymmetry. 
The condition T^j^ > T^^ also requires us to consider the different type of the inflation from 
the previous non-thermal example. 

Since Ki ~ 1 is validated only for the case Tr < Mi, eq. (32) cannot be applied 
to the case (b) in which Tr Mi follows. We can obtain a sufficient amount of the L 
asymmetry by taking Mi large enough to make ei large but keeping sin,^ small enough 
to realize Tr < Mi. Thus, in the case (a) a low reheating temperature like Tr ~ lO^GeV 
can be enough to produce the required B asymmetry by setting Mi and sin,^ suitably. 
To realize such a low reheating temperature, for example, the mass parameters in the 
neutrino sector may be taken as 


Ml = 10® GeV, M 2 = 10^° GeV, M 3 = 10^® GeV, sin^ = 10-^ (33) 

possibility since the reheating temperature is too high to be reconciled with the gravitino problem in this 
case. 

®This oscillation may start during the inflation (F,^ < Mi) or after the inflation (F^ > Mi). 

^^During this oscillation, the flat direction may store the L asymmetry due to the Affleck-Dine mech¬ 
anism as discussed in [23]. However, we do not consider this possibility here. 


16 



Since the effective mass fhi is estimated as fhi ~ 10“®eV, the washout effect is completely 
negligible as expected. We have no gravitino problem in this case since the reheating 
temperature realized by the decay of the Ni condensate is sufficiently low comparable to 
the one given in [21], 

4 Summary 

We have proposed the neutrino mass matrices in the framework of the MSSM extended 
with the three generation right-handed neutrinos. These mass matrices can realize the 
bi-large mixing among the neutrino flavors and explain the neutrino oscillation data. It 
can also saturate the upper bound of the CP asymmetry ei appeared in the leptogenesis. 
Although this model is composed of rather restricted number of parameters, it can make 
the CP asymmetry ei and the effective neutrino mass rhi independent. We have applied 
this model to the thermal and non-thermal leptogenesis and studied the influence of 
this feature on the reheating temperature, which is crucial for the cosmological gravitino 
problem. 

In the thermal leptogenesis our neutrino mass texture seems not to be able to make the 
reheating temperature required from the explanation of the B asymmetry low enough to 
be consistent with the gravitino problem. However, it seems to be able to realize a value 
near the lower bound of the reheating temperature obtained in the thermal leptogenesis 
framework. In the non-thermal case we have found that the low reheating temperature 
consistent with the gravitino problem can be sufficient for the successful leptogenesis. 
Even in that case the parameters in the neutrino mass matrices can be consistent with 
the neutrino oscillation data. 

As is shown in this study, some kinds of neutrino mass texture can be constrained 
by the leptogenesis. It may be worthy to proceed a lot of study based on the concrete 
neutrino model to clarify the relation among the neutrino mass texture, the leptogenesis 
and the gravitino problem. 

This work is partially supported by a Grant-in-Aid for Scientihc Research (C) from 
Japan Society for Promotion of Science (No.14540251). 
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Appendix A 

In this appendix we present an example of the construction for the assumed texture of 
the neutrino mass matrix. We consider the Frogatt-Nielsen type global flavor symmetry 
U(l)^ and an additional discrete symmetry in the lepton sector. The charge assignment 
of the chiral superhelds for the symmetry U(l)^ x Z 2 is assumed as 

w ( 1 , 1 , 1 , 0,1; +), iV 2 (1, 0,1,1,1; +), iVs (1,1, 0,1, 0; +), 

Li (0, 0, -1, -1, 0; +), U (-1, 0, 0, -1, 0; +), L 3 (0, -1, 0, -1, 0; +). (34) 


The Higgs chiral superheld H 2 is neutral for this symmetry. In order to realize the hierar¬ 
chical structure of the mass matrices, we introduce the following several chiral superhelds 
which are singlet for the standard gauge groups: 

01 (- 1 , - 1 , - 1 , 0 , 0 ; -) 02 (- 1 , 0 , - 1 , - 1 , 0 ; -), 03 (- 2 , - 2 , 0 , - 2 , 0 ; +), 

Xi 0) 0) 0) 0;+)) X 2 ( 0 , — 1 , 0 , 0 , 0; -I-), X 3 (0, 0, —1, 0, 0; -I-), (35) 

h (0) 0) 0) 0) +)• 

If we assume that the scalar components of these superhelds get vacuum expectation 
values dehned by 

(36) 


A = M. 

Mpi’ 


(h) 

Mpi’ 


we can obtain both the right-handed Majorana mass matrix and the Dirac mass matrix 
as follows: 


efCV^s eie2CVe3 0 \ / 0 0 0 

AI ~ M 3 eie 2 C^/e 3 ^2.Chz 0 , mo = V 2 (5 (5 , (37) 

^0 0 I ) 0 5 5 y 

where M 3 = Mp\e^ and the order one coefficients are abbreviated. The diherence between 
the cases (a) and (b) should be considered to be explained by these coefficients. 

We can check that these mass matrices can consistently realize the texture assumed 
in the text. Comparing mu in eq. (2) with that in eq. (37) and also using eq. ( 8 ), we hnd 

Applying this result to A4 in eq. (2) and eq. (37), we obtain 


63 ~ 0 . 16 e 2 , sin ^ ~ 



If we take 6 ~ 0.1, ( ~ 0.03 and 62 — 30ei, for example, we find 


Ml ~ 10® GeV, M 2 ~ 10^° GeV, M 3 ~ 10^® GeV, sin^~ 0.03. (40) 


This suggests that (37) can realize the mass matrices assumed in the text. 

We can consider the couplings of the right-handed neutrinos W to the inflaton 0, 
which is assumed to be neutral under the present global symmetry. The lowest order 
superpotential allowed by this symmetry is written as 




(41) 


where the coefficients ci, 2,3 are assumed to be 0(1). Thus the coupling constants \i are 
estimated as 

Ai = ciCiC^, A 2 = C 2 e 2 C^, A 3 = 0363 . (42) 

We hud that these couplings can be consistent with the low value of Tr if ~ lO^^GeV is 
assumed. In this case B 2 3> Bi can be satished since B^ = Ois realized for M 3 ~ 10^®GeV. 


Appendix B 

In the non-thermal leptogenesis the L asymmetry produced through the N 2 decay may 
be washed out in a different way compared with the thermal case. If T^^ < T^^ T^^ is 
satished, N 2 decays into the light particles at the time t 2 immediately after the inhaton 
decays into W- The decay products of N 2 behaves as the radiation and its energy density 
decreases as oc where a is the cosmological scale parameter. On the other hand, Ni 
decays at the time ti after the completion of the inhaton decay. If Ni behaves as a matter 
because of Tr < Mi, its energy density decreases as oc a“®. Then the cosmological 
energy density may be dominated by the Ni energy at least as far as Bi and B 2 are the 
same order. The additional washout can occur in such a case. 

The cosmological energy density p(t 2 ) and the temperature T 2 (t 2 ) of its decay products 
can be expressed as 

= ( 43 ) 
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Taking account of these, we can find that the temperature T 2 of the decay products of N 2 
satishes the relation such as 


(n(t2)y = (LlV'" 

[a(t 2 )) 

From this relation we obtain 

/ -p \ 2/3 

nih) = r.((0 , 


(44) 


(45) 


Now we can estimate the entropy production of the late decay of Ni. Since = 

p(ti)/3Mp[ = Fis satished, the ratio between the entropy density Sb(ti) before the Ni 
decay and the entropy density Sa(ti) after the decay can be written as 

Hih) g.{t2)Tm (m2)y (rf,x 

Saih) g.(ii)rf(ii) tri(«i)/ trj Uw ' 

By using the expression of each decay width, we hnd that k ,2 is written as 


K2 — K 


\/Ml M2 

Tr 


sin^. 


(47) 


where k is the usual thermal washout effect due to the L violating scattering mediated 
by the right-handed neutrinos and so on. 
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